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Ingénierie Scientifique et Visualisation 


[SV 1 DORA Dimension, Optimisation, Réduction et 
Apprentissage pour les problèmes de fluides (Florian De Vuyst) 


ISV2 Ingénierie de systèmes complexes (Christian Saguez) 


ISV3 DIGIPLANTE : Modélisation et simulation pour les 
sciences de l'environnement (Paul-Henry Cournède) 


— ISV4 WYSIWYG : Synthèse d'images - Réalité Virtuelle et 
oimulation (Patrick Callet) 
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Les enjeux de la Recherche en calcul numérique aujourd'hui 


* Physique complexe 

e Aspects pluridisciplinaires 

* Estimation (robuste) de paramètres 

* Assimilation de données 

* Conception / optimisation de systèmes 
* Prise en compte des incertitudes 

e Rendu 


* Temps réel 
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“Physique complexe” - Exemples de problèmes de fluides considérés 


e . Volumes finis pour les fluides compressibles [JCP 2004, ESAIM Proc. 08, E. J. Mech B Fluids 08] 
* Fluides multiphasiques, CEA Saclay [Award ICONE 2004, Nucl. Design 2008] 

e Moulage par injection de thermoplastiques, RNTL, PEP 

*  [nterprénétration de plasmas, CEA DAM {cP 2007] 

e Trafic routier, ACI MD 2003, U. Dauphine, INRETS [Award EGC 2007, Proc. Trans. Res. Board] 
* ` Plasma magnétosphérique, CEA DAM [cPc 2006] 

e Fluides magnétisés, Doshisha U., Kyoto U Phys. Cond Mat., J Mech. Eng . Sci. 2008] 

e Migration/prolifération cellulaire, F. Xavier, Institut Curie, Institut Carnot Supélec 

e ` Modèles réduits pour les fluides, ANR OMD, CSDL [IJNME 2008, Chap. Hermès et Wiley] 

e Biochimie moléculaire dynamique stochastique, Paris 6, ANR CI PITAC [DDM 2008] 

* . Modèles Temps Réel de systèmes (PPF TAMARIS, ENS Cachan) 

e Conception collaborative pour les problèmes à base de fluides (CSDL, Dassault) 

e CO2 transcritique/supercritique (Doshisha U.) 


e Modèles fluides de flux d'informations [Comp. Networks 2008] 
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Coll. France — Japon (+ Chine en projet) 
Conception optimale de systèmes de conversion d'énergie 


CELTE] CELTE] 


Ecoulement de fluide CO2 Ê 
en état transcritique 


m ad (postdoc T. Kuwahara) 
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Devenir des doctorants & postdocs (ISV1) 


* Christophe Audouze Postdoc labo EM2C ECP 


* Fatima Daim CDI IR ENSAM Paris 

* Cédric Enaux CDI Ingénieur CEA Bruyères (ECP 2004) 
* Lucie Fréret CDI IR CNRS Paris 

e Pascal Jaisson PRAG U. Versailles 

* Marc Joliveau Post-doc Montréal 


e Solène Le Bourdiec CDI Ingénieur EDF R&D 
e Takuya Kuwahara  Post-doc ENS Cachan & Osaka U. 
* Michael Ndjinga CDI Ingénieur CEA Saclay (ECP 2003) 
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Un zoom : 


L'initiative 
Complex System Design Lab (CSDL) 


Présentation F. De Vuyst 
Responsable ECP du projet CSDL 


| 
Le projet CSDL 


Mettre en place un environnement collaboratif 
d'aide a la décision pour la conception de systèmes complexes 


Industrie et services 
Grand Groupes PME 
S 


Distene 


Dassault-Aviation 
EADS 


LMS- 
Imagine 


Oxalya 


HPC-Projet 
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Centre de recherche 


Recherche Centre de 
Recherche 


ublique 
ONERA 


Digiteo,/Scilab 


ANSYS ESILV ARMINES 
ESI Group INRIA 


ENS Cachan 


Pilote : Dassault Aviation 


E P=; 


CENTRALE 
Exemple : la conception avant-projet d'avion 
Masse Spécifications avant-projet 
Structure | 
P - | Trainée 
"I. Aérodynamique 


Portance 


Masse de l'engin 


| | Images Dassault Aviation 
Processus de dimensionnement 
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Besoins scientifiques 


e Modèles réduits 


* Propagation d'incertitudes 


* 


systémique 
systémique 


e Modèles de données 


physique 


physique 
physique 


e Requétage complexe 


* Interactions multidiciplinaires 


e Extraction de connaissances 


Visualisation décistonnelle e Visualisation/optimisation “on- 


the-fly” 


Schéma du processus envisagé dans CSDL 
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Projet fédérateur laboratoire + partenaires 


e Marie-Aude Aufaure Modèles et Entrepòts de données 
(chaire Business Objects) 


s Céline Hudelot Recherche de caractéristiques dans les 
grands volumes de résultats, visualisation 


e Erick Herbin Prise en compte des incertitudes 
(+ Supélec + INRIA) Problèmes de probabilité sous jacents 
e Florian De Vuyst Réduction de modèles 


(+ N. Vayatis ENS Cachan) POD et métamodèles pour la visualisation 


Optimisation/Visualisation on-the-fly 
Apprentissage, Scoring / ranking 
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Eclairage : réduction de modèle 


Champ Elément Finis 
Paramètres de conception Spatial (x) 


0 c KR? 


| 


Modele d'incertitudes Y 


u = u(x, 0; Y) 
DU U — u(x, 0: Yo) 


Hyp. : pour une réalisation de Ÿ, on suppose que le champ u est régulier. 
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Modèle de représentation (approximation) 


u(x, 0, Yo) = uo(x) + »- a; (0, Yo) V^ (x). 


A A 


Coefficients de décomposition 
- fonctions de 0 et Yo 


Champ moyen 


Base spatiale 
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Calcul des coefficients 


Si les fonctions V“ sont orthonormées 
dans L“(Q), alors 


a; (8, Yo) — (nt. 0, Yo) — uo(-), V*(-)) r2tgy: 


/ 


Processus stochastique a déterminer numériquement 
(polynómes de chaos, ...) 
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Choix possible de base réduite spatiale 


Les fonctions V^ sont choisies comme modes POD. 
Méthodologie : 


1. Plan d'expérience numérique explorant O et Y o 
2. Base de données de solutions 
3. ACP sur la base de données (choix d'un produit scalaire) 


Approche non intrusive 
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Discussion 
dans le cas déterministe 


u(x, 0) & uo(x) + X p=. a1(0) hier), 


a 


"a 


P 


Etude de cas - Solutions de Navier-Stokes dépendant 


de onze paramètres 


[Audouze, De Vuyst 2008] 
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> Total de 11 paramètres CAO 


Chaque paramètre varie dans un intervalle de conception 
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Plan d'expérience numérique — Exemples de solutions EF 
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Plusieurs types de solutions 
Comité d'Evaluation MAS - 4 décembre 2008 


18 


x 
FF = me E 
Ki 1 
d É 
S J 
a 
po] 
- L 
mr > L 
ek # 
= 
" ML 
k. 
J L 
T nil h. 
L 


oolutions discrètes EF paramétrées 


Espace paramétrique : 
u” = u” (x, 0) 


0 c |0, 11” 
/ \ 
Variable Vecteur ` 
Q : comment construire un estimateur “rapide” de 


9 u"(.,0)? 


( naz Comité d'Evaluation MAS - 4 décembre 2008 19 


LELE 
# 


+. 


EEK 


eg 

ep e d 

- L 

-* ^, 

ans 
Ze 
D 
Oh 
` 
` 
` 
Na 
paas 


` 
-"—— => 


-- 
=... 
d ës ee inn mea a 


“e 


an 


— XN N X XX X 
- 


IM 


e — ^ 
CA =` ` 
LELE EEE E a 
ee mo-a 
"T2227 mm SN kA 
TIP FD ma ANA 
LÉ x NN AÁ 
#esgezse-sssii 


èn 


is 
s- 
Spaa 


9 


` 
. 


Bou ott mS 
` 
x x x X x x x X man 


^ 
x as K «œ --— 


- 
- 
- 


ns -# 
"a 


rèn AAA 


` 
997 ò ta => «< vann 


4044 
EZ 


ees ER 8 8 o wm 


“sssa 


e LA EPE ECE e 
Sare r PC y statut 
papan 

sa 44 4 44 a 507 
«ennek KX X XX 
ressvsdsèar minn 


4442229» 999 ée ees «œ sm 


` 


eanna XM PA PA 
X XXX XNXX a m r e e ée e zm 
XX A TU m ne rr e a P d rus C 
#' 
---. 


H 
H 
` 
` 
` 
` 
` 
H 
T 


À 
ar pp sn ees ee ee è eler 


+ 
D 
, 
D 
a 
» 
D 
` 
` 
` 
a 
D 
1 
` 
` 
` 


H 
`` 
“a 
«x 
>x 
=. 
e x 
*. 
vé 

` 
L. 
` 


D 
D 
U 


rer e opp pm ana an a 


sann 


—o-—— È mem e ee —— nn 


X XXX N ON Ten e emen an sas omm as =W" ` 


yo 
e 99999 = ee eege e eege Ayi ponn. 


- 


a sèn Ta 
sa 


LEA 
CI att. 


of ont * "^" ^^ 


x XX X N N mm 
X XXX Sr em are er an 
TL 
\ 4044064060606060 6 aae =œ K Y N X X K 
dottisseceserrepam- atat 
t 404 444 tt 


a seo AA AJI ENN “au 
et" # 


res 
LED 
ER 


0292209609 20 


444 ## AA 


` 
D 
H 
` 
` 
` 
` 
` 
` 
` 
` 


UREE E KAN a 
rt eg 
a 

sz seleeeg 


1IIIIIIIIEIIE d fraen 
D 


KEK CS NC NTC UL H da 


17:33 DN 
"9 89 


""xxazx"PPPFFE'x.H 


x + 
LK ITR T 
` 
` 
a a sn an 


an T APA H R 
` 


rr 
XX rat na Y K X Y 


XXX N NNR AAA X K KKR E kd 
d ee wies eê 6 


s. 
sisan Mi 
X X X X Xx 
XXX X XX XxX X Xx 
Sx s N X YX xY x X 
=s XNXX > + 


wes SAN LKKL KAN 


sa 
XXX RMX RY mr 


ta 

isa. 

veso van n nn 25 P488 
x x n 

FÒ ce v 


H 


«ee X XXX ` 
sss XXX YX x 


e wë TT ET mg ep e d d 


D er em 2" mp 


EA Ba H N 


` 
` 
» 5 
“un 
a` 
+ 
ra 


L 

` ss... 
bn ans 96 AD "^ === $ à # 4 
`... 


- 
` 
' 
* 
D 
D 
a 
a 

e D 

-- - 
- 
D 
a 
D 
H 
D 
D 
` 
` 
D 
e 
- 


Aan 
CE S ` 
``. 
Vous 


mr, 
“ts 


POM #1 


89 
382 
339 


œ 
© 
© 
av 
6, 
Da 
fe 
= 
d 
O 
`Ü 
O 
<T 


V) 
< 
2 
c 
Kë 
Œ 
= 
Œ 
> 
LII 
O 
`O 
E 
O 
O 


[ON 


A 155 
À" 145 
85 


de la 


T 


ecroissance 


en vitesse 
A 
D 

Oe corrélat 


3 


95 
854 
24 


. 
` 


lères v.p 
matrice 


Modes POD 


23708 
A? 3667 

735. 

505 


A? 
M 


AM 
Bon taux de d 


a nn 


Ü Í i sijè 


Comparaison solution E.F. / approximation POD 


GA 


renn 


— 


N 


“SNN 
SSN NYK 
e lab. A 


/ 
1 
/ 
U 


m 
P ortum mn 


MW ét am 
f ë P wes 
VIL dè m 


E 


ti P 9 ame 
Ka f am 


mL sa su DE trann Sa 
MT gi R — 


A 
/ 
/ 
li 
Í 
| 


N mm N a Ta Dam non t 


Reconstruction avec K=5 


(5 modes POD) ! 
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Apprentissage des coefficients 


1. Generer un nuage de N points 0" dans 0, 1}? (LHS, Sobol, etc. ); 
2. Calculer N solutions discrètes EF yu!" = u? ^^. 
3. Calculer a} = a, (0°) = (ach, y^^: 


4. Interpoler/approcher les coefficients a; dans le domaine parametrique. 


[Audouze, De Vuyst IJNME 2008] 


> Nécessite des interpolateurs/approximateurs haute dimension 


Plusieurs candidats : 


Polynómes 

Fonctions à base radiales (RBF) [Wendland 2006] 
Approches à krigeage 

Moindres carrés mobiles (diffuse approximation) 
Approximation particulaire type SPH, etc... 
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Expérimentation numérique 


Système non linéaire 1D de Keller-Segel 
(modélisant la formation de motifs sur la peau 


de certains étres vivants) you Termes non linéaires 


Ou — Au + aV - (uVc) = Eu H — u) 


d 
Oe — DAC = — c| 
i 1 ( ) 


Structure riche des solutions dans l'espace paramétrique. 
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Solutions 1D stationnaire avec CL périodique 
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Paramètre y variable seulement dans cet essai. 
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4 Décomposition POD 


D: 


Correlation matrix Eigenvalues 


OO NN 6H A Aa 


Spectre de la matrice 
de corrélation 
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Les huit premiers modes POD 
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CENTRALE Comparaison POD-RBF ROM / Eléments finis 


(en utilisant K=6) 


EE z 
9.0.3.40.0-8 .Œ S ©.2 4 .@ 81.0 7 H 7.055 A 2 


103404 0.2 © 3 0 © 4 0 Ke 005.098 0.2 


© Mo LH GO n CAMION) CLRH 


Légère perte de précision 

au voisinage de la bifurcation 
mais restitution globale 
satisfaisante 
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Construction pratique d'une 
partition de l'unité : 


T 


o^ 


ò" (9) = d 


Noyau régulier positif 
a support compact 


- w” (0 | 
Alors w™ (0) = ——— —— ales propriétés attendues : 
^d 
` W (9) m M = 
TE W = 0, ) m W = Ï. 
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Algorithme 


1. Construire une base locale POD relative 
a la fonction de partition n? m : 


(N | K m : 
+ A 
min — A wm(0) Ki TAN (ab pwe pin 
BEE 2 = | i 
(WE L^ )—6L,, 1SKSLE KA 
li Moindres carrés a poids 
K m 
= — 0 bh N : h 21 k. h k. h 
Reconstruction locale: U ^ = » (u (8 ). Vo ) Wo 
i k=1 
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e 
Reconstruction globale 
e A partir de 


» wm (8) = 1 


on a l'approximation globale : 


M ja Gon 
au h kh\ atik, 
m > Wm (8) (u (9) Vi ) Wi” 
m=] k=1 
de la forme Kg Approximateurs KBF 
r = Y CL 
m= k=1 
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Un élément d'analyse numérique 


(cas d'un problème d'évolution) 
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l 
Approximation POD 


x 
FF = me E 
"E i 
"d P J 
S J 
a 
po] 
- L 
E a Ï 
„4 # 
s= 
" ML 
k. 
Tu L 
T nil h. 
L 


Most of the information 
contained here 


Soit 
e > U ` 
un seuil d'erreur. Si k 
K; / S` à" < 2» (Ky = Kÿ(e)), Eigenvalue 
EN residual 


alors on a l'estimation part 


ol IK: fo < SÉ 
7= || 
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Ex éq. de la chaleur 
“Réponse” de la solution a la donnée 
initiale et au second membre 


TI p 
Te | = oi Wf, 
Fi a F 
Be | 
L =| 1 L 
= BM 
_ | Tœ WA 


= (u, Vine F (u,v)1,0,0 = (f, v)o,o g(t) ` for all v € H;(Q), te, TI) 


u(,t=O)=u inn. (2) 


SY ZIP... e 1 (Q)) 
$54 = 10 w € L (Q)V. 


uo — 


Kj Ko 
K 
H" f => (f, V5)oo V5, f ELO), IIg^ u? => (u, Vo)oa Yo, u? € LAA). 
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S 
CENTRALE Dual basis construction 


— (p4, v)o,o T (47. v)1,0,0 — (Pi, vog g(t) for all v € H5(Q),t = (0, T), 


dt 
QL, t= 0)=0 inQ 


d 
(po: Ua. 0 an (26, v)1,0,0 — U for all v € H (Q), l E (0, T), 


po.,t=0)=VW inn 
The solution ü of the approximate heat equation problem 
d ,- : 
mu v)o,q + (ü,v)1,0,0 = (II; 7 f,v)oa gt) forall oe HA, (A), t € (0, T), 
a, ,t=0)=I,u inn 


writes in closed form 
Kf Ko 


ü(-, t) — » (f. Yon ert) T N (u, V 5)o,o Pol, t) Vt € 0, T'). 


Ç ad Comité d'Evaluation MAS - 4 décembre 2008 33 


l " 
L 


Estimation d'erreur a priori 


En utilisant le lemme de Gronwall, 


Ye tegere dont 


with A(t) = Leen rés € z) g(s)|ds 
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Deuxième zoom : 


Projet DIGIPLANTE 


Resp. P. H. Cournède 


